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Abstract. Backscattered ultraviolet (BUV) instruments designed for measuring stratospheric 
ozone profdes have proven to be robust tools for observing polar mesospheric clouds (PMCs). 
These measurements are available for more than 30 years, and have been used to demonstrate the 
existence of long-term variations in PMC occurrence frequency and brightness. The Ozone 
Monitoring Instrument (OMI) on the EOS Aura satellite provides new and improved capabilities 
for PMC characterization. OMI uses smaller pixels than previous BUV instruments, which 
increases its ability to identify PMCs and discern more spatial structure, and its wide cross-track 
viewing swath provides full polar coverage up to 90° latitude every day in both hemispheres. 

This cross-track coverage allows the evolution of PMC regions to be followed over several 
consecutive orbits. Localized PMC variations determined from OMI measurements are 
consistent with coincident SBUV/2 measurements. Nine seasons of PMC observations from 
OMI are now available, and clearly demonstrate the advantages of these measurements for PMC 
analysis. 


1. Introduction. 

Satellite measurements of backscattered ultraviolet (BUV) radiation have been used to 
monitor stratospheric profile ozone and total column ozone since the late 1960s [ Heath et al., 
1973]. Regular ozone observations began in November 1978 with the Solar Backscatter 
Ultraviolet (SBUV) and Total Ozone Mapping Spectrometer (TOMS) instruments on the 
Nimbus-7 satellite [Heath et al., 1975]. Additional TOMS instruments were flown on the 
Meteor-3 [Gleason et al., 1993], ADEOS [ Krueger et al., 1998], and Earth Probe [McPeters et 
al., 1998] satellites through 2004. A second-generation version of the SBUV instrument, 
designated SBUV/2, was developed by NOAA for operational use. Seven SBUV/2 instruments 
have been flown since 1985 [e.g. Frederick et al., 1986; Hilsenrath et al., 1995; DeLand et al., 
2001], and four SBUV/2 instruments are currently operating. 

Polar mesospheric clouds (PMCs) are tenuous, narrow layers of ice crystals that form at 
80-85 km altitude only during a well-defined season in each hemisphere (approximately 20-30 
days before summer solstice to 60-70 days after solstice) [e.g. Thomas, 1995]. PMCs are 
generally observed at latitudes poleward of -50°, although observations down to 40°N do occur 
[e.g. Taylor et al., 2002], Thomas et al. [1991] demonstrated that Nimbus-7 SBUV 



measurements of backscattered albedo, at wavelengths between 252 nm and 293 nm, could be 
analyzed to identify the presence and brightness of PMCs in the upper atmosphere. Since these 
measurements are taken during the normal operating mode of the SBUV instrument, PMC 
information can be derived from data obtained by all subsequent SBUV/2 instruments as well. 
PMC results from combined SBUV and SBUV/2 data sets have been presented in numerous 
papers [e.g. DeLand et al., 2003; Thomas et al., 2005; DeLand et al., 2007; Shettle et at., 
2009], These results have demonstrated the existence of two significant long-term variations in 
PMC occurrence frequency and brightness: Increasing trends at most latitudes in both 

hemispheres, and an anti-correlation with solar activity over solar cycle time scales. However, 
our ability to evaluate SBUV/2 measurements for PMC properties and morphology is 
constrained by the nadir viewing geometry and large field of view (nominal size = 1 1 .3° x 1 1.3° 
= 170 km x 170 km at the Earth’s surface). The SBUV/2 instrument samples wavelengths 
sequentially with a 2-second interval during normal operations, so that the field of view moves 
an additional 70 km along the orbit track during the measurement of the five wavelengths used in 
PMC analysis. The complex small-scale PMC structures shown by data from the Cloud Imaging 
and Particle Size (CIPS) instrument on the Aeronomy of Ice in the Mesosphere (AIM) satellite 
[Rusch et al., 2009] clearly shows that SBUV/2 measurements typically observe both PMC and 
non-PMC regions in every sample. This averaging process thus reduces the ability of SBUV/2 
measurements to detect faint PMCs, as well as narrow and isolated bands of bright PMCs. 

The Ozone Monitoring Instrument (OMI) is a hyperspectral instrument flying on the EOS Aura 
satellite to continue the stratospheric profile ozone and total column ozone record [Levelt et at., 
2006a]. Aura was launched on July 15, 2004, and regular operations began in December 2004. 
OMI measurements have a smaller pixel size and wider geographic coverage than SBUV/2 
measurements, offering significant advantages for PMC detection and analysis. Approximately 
30-40 OMI nadir view pixels fall within a single SBUV/2 footprint, improving the ability to 
identify small-scale PMC structure. OMI observes most locations at typical PMC latitudes 
(>55°) on two or more consecutive orbits each day, and its coverage extends to 90° latitude in 
both hemispheres. This paper presents an overview of the OMI instrument and the adaptation of 
the SBUV/2 PMC detection algorithm for use with OMI data. 

2. OMI Instrument and Data Description. 

Detailed descriptions of the OMI instrument design and operation are provided by Levelt 
et al. [2006b] and Dobber et al. [2006], OMI is a single monochromator that covers the 
wavelength range 264-504 nm, separated into three spectral ranges with overlapping coverage. 
The measurements used for PMC analysis are taken from the UV1 channel, which covers the 
region between 264-311 nm with 0.63 nm spectral resolution. OMI views the Earth in the nadir 
direction using a wide field-of-view (FOV) telescope with a total cross-track coverage of -1 14° 
(Figure 1). The UV1 data are binned into 30 rows with equal angular extent, yielding pixel 
cross-track widths of 48 km at nadir (rows #15-16), increasing to -200 km at the outer edges 
(rows #1 and #30). OMI co-adds five consecutive measurements with an exposure time of 0.4 
seconds per measurement, which gives a nominal pixel size of 13 km (along-track) x 48 km 
(cross-track) in the nadir rows, increasing to -25 km x 200 km at the outer edges of the viewing 
range (Figure 2). We do not use rows #1 and #30 in this analysis due to problems with signal-to- 
noise performance. Over the course of OMI’s 14 daily orbits, the total cross-track swath of 
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approximately 2600 km provides nearly full global coverage at the Equator, with substantial 
overlap between consecutive orbits at PMC latitudes. 

Instrument throughput changes in the UV1 spectral region have been analyzed using 
polar ice cap radiances [Jciross and Warner, 2008], The OMI calibration has been shown to be 
stable to better than 0.5% [G. Jaross, private communication ]. Since OMI is a single 
monochromator, the radiance measurements at wavelengths less than 300 nm can have 
significant stray light contributions from longer wavelengths because of the rapid spectral 
variation in backscattered radiance [e.g. DeLand et al., 2001]. The operational OMI data employ 
a stray light correction algorithm developed from both prelaunch calibration tests and inflight 
measurements using dedicated rows of the CCD. Dobber et al. [2008] describes the use of both 
spatial and spectral stray light corrections in the creation of OMI Level lb data. These 
corrections are estimated to reduce the stray light contamination for UV 1 measurements to less 
than 5%. 

We use OMI Level lb radiance data from the Collection 3 processing in this analysis. 
These data include quality flags for a variety of parameters at the pixel, swath, daily, and global 
levels of instrument data [van den Oord et al., 2006]. A specific flag of interest is the random 
telegraph signal (RTS) flag, which indicates raw pixels whose dark current signal tends to shift 
between multiple stable levels as a result of proton impact damage [Dobber et al., 2006], 
Approximately 30-50 raw pixels are averaged into each albedo bin used for PMC analysis. 
Detailed examination does not currently show an impact on OMI PMC retrievals when RTS 
flags are present in the lb data. However, we do save these flag values in our PMC data sets, 
and we will continue to track the behavior of the retrievals. 

OMI Level 2 total ozone data began to show row-dependent anomalies in June 2007. 
These anomalies were soon identified in Level lb UV1 radiance data as well, and have continued 
to evolve with time. The detailed behavior of the row anomaly is complex, but the general effect 
is to enhance observed radiance signals in affected rows in the Northern Hemisphere PMC 
region (latitude > 40°) by factors up to 5-1 Ox, and to decrease observed Southern Hemisphere 
signals in affected rows by factors of 0.6-0. 8. The magnitude of the observed effect also varies 
with solar zenith angle (SZA). Detailed discussions of the row anomaly behavior can be found 
at http://www.knmi.nl/omi/research/product/rowanomalv-background.php . At this time, no 
correction process is available for these effects in UV1 data. Based on inspection of PMC 
analysis results, we have chosen to omit any row from our data processing if its background 
deviation from the SZA-dependent variation of the unaffected rows is greater than 10%. This 
represents rows #16-23 for the NH 2008 season, and rows #14-23 for the SH 2008-2009 season. 
Initial NH 2009 results show an expansion of the contaminated region in UV1 data to cover rows 
#7-26. We will continue to track the status of instrument throughput changes and row anomaly 
behavior. 

3. OMI PMC Detection Algorithm. 

The measurement technique and spectral coverage of OMI allow us to adapt the SBUV 
PMC detection algorithm [Thomas et al., 1991 ; DeLand et al., 2003, 2007] to analyze OMI data. 
This algorithm looks for pixels with albedo enhancements relative to the Rayleigh scattering 
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background using data at the five shortest SBUV wavelengths (252-292 nra), and uses a 
combination of four tests to distinguish pixels containing PMCs. Since OMI spectra have a 
slightly different wavelength scale for each swath, we interpolate all spectra to a regular grid 
with AX = 0.5 nra for consistent processing. OMI spectral coverage does not extend as low as 
SBUV/2 instruments, so we adopted the revised wavelength set 267.0, 275.0, 283.5, 287.5, 292.5 
nm for OMI use. We plan to conduct further studies to investigate possible benefits of using 
different (or expanded) wavelength sets. We found that averaging albedo values for three 
consecutive 0.5 nm bins at each PMC analysis wavelength significantly improved the noise in 
the albedo data. 

The Rayleigh scattering background atmospheric albedo observed by OMI at the 
wavelengths used for PMC analysis varies by almost an order of magnitude with changing solar 
zenith angle, and the off-nadir measurements typically show increased albedo values relative to 
nadir measurements by up to a factor of two due to the increased viewing path through the 
atmosphere. In order to allow processing of all OMI data with a single set of algorithm 
parameters, we scale the data for each row to be consistent with nadir-viewing measurements 
using geometric terms which are independent of the variable atmospheric properties. The 
following discussion presents the derivation of this scaling procedure. 

The background atmospheric albedo, /Cay, due to molecular Rayleigh scattering in the 
presence of ozone absorption is given by Bailey et al. [20091, who generalized an earlier result of 
McPeters [1980]: 

A Ray = PRay(©)r(0+ 1 )pMz 0 ) / { Pv[( 1/pv) + Chi gsz A H oz )] a a rT C( Zo) 0 } [ 1 ] 

P] R ay (©) is the Rayleigh scattering phase function for a scattering angle ©, T is the gamma 
function, o is the ratio of the ozone scale height (H oz ) to the atmospheric scale height (assumed 
to be constant for the altitudes of interest), p is the cross section for Rayleigh scattering, N is the 
column density of atmospheric molecules above a reference altitude Zo, pv is the cosine of the 
view zenith angle, Psza is the cosine of the solar zenith angle, Ch is the Chapman function 
[Chapman, 1931], a is the cross section for pure absorption by ozone, and C is the column 
density of ozone molecules above that reference altitude. For solar zenith angles less than 70° 
the Chapman function reduces to 1/pszA- 

Equation [1] for the atmospheric albedo can be separated into two factors, /(Zo) which 
just depends on the atmospheric properties and G( a v- Psza) which depends primarily on the 
viewing/sun geometry. These factors are defined as: 

F(z 0 ) = Ho+DP^zo) / [a CT C(z 0 ) a ] [2] 

and 

G(bv> PsZa) = PRay(©) / {hv[(l/bv) +C/l(psZA,H oz )] C '} [3] 

It will be noted that the geometric factor, G, still has some dependence on the atmospheric 
properties through the ozone scale height, H oz , and exponent a, the ratio of the ozone scale 
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height to the atmospheric scale height. However, for the wavelengths of interest (267 to 292.5 
nm), ozone absorption limits the backscattered radiance viewed by the satellite to altitudes 
between approximately 40-60 km. At these altitudes, the typical values of cr vary from 0.5 to 
0.65, and the use of a median value between 0.55 and 0.58 in Equation [3] reproduces the albedo 
variation with solar zenith angle to 1 0-20%. These variations are due primarily to uncertainties 
in cr, which is taken to be 0.58, rather than uncertainties in H oz affecting the value of the 
Chapman function. Even a factor of two change in the ozone scale height changes the geometric 
factor, G, by as much as 10% only for our most extreme solar zenith angle of 88° Dividing the 
observed atmospheric albedo A 0 bs by the geometric factor G given in Equation [3] then removes 
most of the variation of the atmospheric albedo due to changes in the scattering geometry. 


Ascaled ~ A 0 bs / G(,Uv, PsZa) [4] 

This makes it easier to approximate and remove the background atmospheric radiance from the 
measurements with the clouds. Using this procedure, we find cross-track agreement in absolute 
albedo between rows of OM1 data to better than 10%, and gives adjusted albedo values for each 
row in a format and magnitude consistent with SBUV/2 data. The agreement in albedo between 
rows reduces the uncertainty of the background albedo fit to 1-2 %. 

For SBUV/2 data, measurements from all orbits for a single day are combined before 
calculating the background fit for PMC analysis. The higher along-track sampling rate of OMI 
compared to SBUV/2 (13 km vs. ~220 km) provides sufficient measurement density to process 
each OMI row separately, and furthermore allows each OMI row to be analyzed individually 
within a given orbit. This allows the background fit to the adjusted albedo to be more responsive 
to any longitudinal variations in stratospheric ozone. We used out-of-season measurements to 
represent typical atmospheric variability, and calculated the albedo standard deviation as a 
function of latitude at each wavelength. A quadratic fit to these data provides a threshold 
function that is used to screen out anomalous samples, and is also used to define the residual 
albedo threshold at 267 nm for PMC detection (Figure 3). Note that the relationship between 
latitude and solar zenith angle varies during the PMC season. We found that scaling this 
threshold function up by a factor of 1.6 was required to decrease PMC occurrence frequencies 
below 1% at the start and end of our nominal season (30 days before summer solstice to 70 days 
after solstice) [DeLand et al., 2003], as expected based on results from more sensitive satellite 
instruments such as SNOE [Bailey et al., 2005], The combination of higher spatial resolution 
and greater geographic coverage means that OMI detects many more PMCs than an SBUV/2 
instrument, which greatly improves our statistical basis for examining PMC behavior. 

4. Sample Results. 

Figure 4 shows an example of the OMI PMC albedo values retrieved for orbit 15878 on 
July 10, 2007 (day 191). Coherent structure is visible between 70°-80° latitude in numerous 
places, including a large bright region at 68°-73°N and 75°-95°W. Figure 5 shows the residual 
albedo for row #17 only, which extends from approximately 40°N, 50° W to 65°N, 170°E in 
Figure 4. Structure at the sub- 100 km scale is evident at SZA = 70°-80°. Note that some pixels 
with residual albedo values exceeding the threshold shown in Figure 3 are not classified as 
PMCs based on the results from additional detection tests, as described in DeLand et al. [2003, 
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2007], In some locations of this orbit, retrieved cloud albedo values can exceed 50 x 1 0 6 sr 1 . 
This is brighter than almost any PMC albedo value observed by an SBUV/2 instrument, and is a 
natural consequence of the fact that it is easier to fill a 750 km 2 OMI pixel with PMCs than a 
22,500 km 2 SBUV/2 pixel. OMI off-nadir measurements can have significantly different 
scattering angles (up to 20-30 degrees higher or lower) than the coincident sub-satellite 
measurements, so that ice phase function variations may affect the observed PMC brightness 
[Bailey et al., 2009], Some of the residual albedo values are negative because of the remaining 
uncertainties in our fit of the background, which is subtracted from measurements, as may be 
noted in Figure 5. 

The orbits of the Aura (with OMI) and NOAA-18 (with SBUV/2) satellites have very 
similar Equator-crossing times at different altitudes (and therefore have different orbital periods), 
so that Aura passes beneath NOAA-18 every 31 orbits, providing good spatial and temporal (A t 
<15 minutes) coincidences [ Shettle et al., SBUV observations of polar mesospheric clouds 
compared with MLS temperatures and water vapor measurement, in preparation]. For these 
situations, we find that selecting OMI pixels whose center lies within 75 km of the center of the 
SBUV/2 field of view (FOV) gives approximately 25-28 OMI pixels within each SBUV/2 FOV. 
Figure 6 shows a comparison of the NOAA-18 SBUV/2 residual albedo values at 252 nm for 28 
SBU V/2 events on July 3, 2007 with the average of all coincident OMI residual albedo values at 
267 nm, inside the FOV for each SBUV/2 measurement. Squares indicate the samples identified 
as PMCs by the SBUV/2 algorithm. The correspondence is reasonably linear despite the 
relatively small number of data points, with a correlation coefficient of 0.87. The majority of the 
PMC detections have higher albedo values for OMI than for SBUV/2, consistent with our 
previous statement that OMI is capable of identifying brighter clouds. 

The large geographic coverage of OMI measurements makes it possible to track bright 
PMC regions over successive orbits. Figure 7 shows the retrieved PMC albedo for six 
consecutive orbits on 2007 day 191 (orbits 15873-15878), focusing on the region 60°-90°N. A 
bright region is visible to the west of Greenland at ~70°N on orbit 15874, and moves to the west 
on subsequent orbits. The brightness and internal structure of the region also evolves during this 
time. NOAA-17 and NOAA-18 SBUV/2 observations on this date also show the presence of a 
bright region of PMCs near Greenland. The orbit-by-orbit OMI observations also indicate the 
presence of extended voids at high latitudes ( e.g . 80°-90°N and 1 10°-140°W in orbit 15878), as 
first noted by Rusch et al. [2009], 

Figure 8 shows the seasonal variation of PMC occurrence frequency observed by OMI in 
5° zonal average latitude bands for the NH 2007 season. Relatively few PMCs are seen below 
~60°N. Daily average occurrence frequency values of 40-50% are commonly seen above 80°N, 
again higher than corresponding values for SBUV/2 instruments. Figure 9 shows the 
corresponding seasonal frequency plot for the SH 2007-2008 season. PMC frequency values are 
lower than NH values for latitudes less than 70°, but become significantly higher above 80° 
latitude. It should be noted that typical scattering angle values at 85°S range between 80°-90°, 
whereas scattering angle values at 85°N are approximately 110°- 120°. The ice phase function 
tends to be smaller at larger scattering angles, although the magnitude of the difference depends 
on details of the corresponding particle size distribution. Thus, increased PMC albedo values 
due to ice phase function effects may lead to higher observed occurrence frequency values in the 
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Southern Hemisphere high latitude data. A complete interpretation of the PMC latitude 
dependence as observed by OMI also requires an understanding of the PMC local time variation, 
which is discussed further in DeLand et al. [2010]. 

5. Conclusions. 

The OMI instrument represents an important advancement in satellite observations of 
PMCs. While its basic design provides continuity with the multi-decade SBUV series of 
instalments, the optical design and detector characteristics of OMI greatly improve its usefulness 
for PMC analysis. The use of a 2-D CCD detector produces smaller spatial pixels with 
continuous UV spectra for every pixel. The wide cross-track measurement swath provides daily 
observations up to 90° latitude in each hemisphere, and gives multiple coincident samples every 
day at most PMC latitudes. DeLand et al. [2010] describes how we utilize this expanded 
temporal sampling to characterize PMC local time variations. In the future, we intend to use the 
spectral coverage of the OMI measurements to derive Angstrom exponent values and thus 
investigate PMC particle size behavior. 
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Figure 2. OMI pixel size and row location for UV1 wavelengths relative to orbit track. Note 
that the cross-track and along-track plot scales are different. 
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OMI Residual Albedo for 2007 / 191 : Orbit = 15878, Row = 17 
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Figure 5. OMI albedo residual values [observed - background] at 267 nm for all samples from 
row #17 of orbit 15878 on July 10, 2007, plotted as a function of solar zenith angle. Squares 
represent pixels where a PMC was detected, and crosses represent no PMC detection. The 
corresponding latitudes for regularly spaced locations are identified with tick marks at the top of 
the plot. 
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Average OMI Albedo vs. SBUV Residual Albedo 



Figure 6. Residual albedo values for coincident NOAA-18 8BUV/2 and OMI measurements on 
July 3, 2007. Coincidence criteria are defined in the text. Squares identify the samples where a 
PMC was detected in the NOAA-18 sample. The dashed line indicates a one-to-one relationship 
in residual albedo. 
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Orbit = 15373 


OMI Data for 2007 / 191 : Consecutive Orbits 

Orbit = \ 5874 
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Figure 7. PMC albedo values at 267 nm observed by OMI for 6 consecutive orbits on 2007 day 
191. Data are shown for 60°-90°N latitude. Note the evolution of the bright region that is first 
observed at 70°N, 45°E in orbit 1 5873. 







